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A persistent current qubit has two quantum states with 
opposite currents flowing in a superconducting loop. Their 
magnetic field couple to nuclear spins. The qubit state is not 
only perturbed by the spins but it also gets entangled with 
the spins' state on a very short timescale. However, when 
the same spins are exposed to a strong but less than critical 
external magnetic field, then the qubit field is just a small 
perturbation on top of the external field and the entangle- 
ment with each spin is negligible. For a qubit which is more 
microscopic than certain threshold this partial entanglement 
results in negligible partial decoherence. 



I. INTRODUCTION 

The idea to perform computation with the help of 
quantum mechanics dates back to Feynman [jlj . Research 
in this direction got a considerable acceleration after the 
realisation (5) that quantum computers could solve in a 
polynomial time certain problems which require nonpoly- 
nomial time on classical computers. Two main obstacles 
to a universal quantum computer are decoherence and 
scalability. Decoherence, if sufficiently weak, can be dealt 
with by quantum error correcting codes, see Ref. @ for an 
experimental implementation. Scalability requires that a 
given technology can be upscaled to a quantum coherent 
circuit with hundreds or thousands of quantum gates, 
see Ref. Q for reviews on various proposals of scalable 
implementations. At the present stage solid state imple- 
mentations have most to offer from the point of view of 
scalability but at the same time they are most likely to 
suffer from decoherence. The aim of this note is to re- 
duce decoherence in the persistent current qubit (PCQ) 
proposal g. 

Two remarkable experiments Q took place last year 
where a superconducting loop broken by Josephson junc- 
tions was forced into states which are quantum superpo- 
sitions of two states with opposite macroscopic persistent 
currents: | f) and | J.). The two considered combinations 
|±) = I T) ± I J.) were eigenstates separated by a gap of 
OAK so it is not surprising that each of them does not 
suffer from decoherence at the temperature of 50mK. A 
more general superposition a\+) + /3\— ) would quickly 
decohere into a mixed state |a| 2 |+)(+| + |/3| 2 |-)(-|- To 



use a persistent current qubit for quantum computation 
any superposition must be free from decoherence. A de- 
sign of a more microscopic scalable PCQ was proposed 
in Ref. ||. These authors went on to estimate decoher- 
ence times from different sources relevant for their design 
and found out that the most dangerous are nuclear 
spins which couple to the magnetic field induced by per- 
sistent currents. Similar conclusion is reached in Ref. 

In the calculations of Ref. the spins are treated 
as random static background fields as may be justified 
by the long relaxation time of nuclear spins which is of 
the order of minutes. Such a static background can be 
measured at the beginning of the quantum computation 
and compensated for by an adjustable counterterm in the 
qubit Hamiltonian. In this note we re-address this prob- 
lem and find out that spins cannot be treated as a mere 
random background because the qubit state quickly gets 
entangled with the spins state. This observation is in 
agreement with Ref. B. What is more, because of the 
large number of involved spins, even for the long relax- 
ation time the spins background performs a random walk 
which leads to dephasing in the qubit state. The former 
more acute entanglement problem can be treated by ex- 
posing the nuclear spins to a strong external magnetic 
field B parallel to the plane of the superconducting loop 
and penetrating through its thin layer, see Fig.l. 




FIG. 1. The idea of how the external magnetic field B 
can be introduced into a superconducting loop. Only one 
Josephson junction (JJ) is marked. 

B has to be less than the critical field in a type I su- 
perconducting film. If the N nuclear spins are exposed 
to a magnetic field b from the persistent current which is 
much weaker than B, as can be characterized by a quality 
factor 
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then the entanglement between the qubit state and the 
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spins state is negligible and leads to negligible partial de- 
coherence. The means which treat entanglement inciden- 
tally also reduce the latter random walk problem, which 
can be furter eliminated by increasing the spin relaxation 
time T r which depends exponentially on temperature. 

II. SINGLE SPIN ENVIRONMENT 

We begin with an elementary example where the qubit 
is coupled to just one environmental 1/2 spin. For the 
purpose of this paper the spin-1/2 is a sufficient repre- 
sentation of the nuclear spin-9/2 of niobium. The spin is 
coupled to the external magnetic field B; its Hamiltonian 
is 

H E = Baf . (2) 

The current of the qubit induces a magnetic field which at 
the location of the environmental spin has strength b. For 
the time being we assume that the field is perpendicular 
to the external field in, say, x-direction. The Hamiltonian 
of interaction between the qubit and the spin is 

ffq-E = ftowtrf . (3) 

The external magnetic field B, which is in plane of the 
superconducting loop, does not couple to the qubit state. 
For the sake of simplicity we set the qubit Hamiltonian 
Hq = 0, what does not change the general conclusion of 
this note. The initial state of Q and E is a product 

|V(0))(a|T>+/?|4»UE>. (4) 

The states | 1(E)), I T(E)) are ±1 eigenstates oi a z . Un- 
der the total Hamiltonian 7?q_e + He the initial \ip{0)} 
evolves into 

Mt)>=a| T) e - lt(+ ^ E+S ^ ) llE) + 

P\ i)e-^-^+^)\ iE ) . (5) 

At the same time the initial pure reduced density matrix 
of the qubit 

3?) <6) 

evolves into 

M - ( a * a 0(t)a(3* \ 
p[t) ~ \0(t)a*[3 (3*f3 ) ' {t) 

where the overlap 0(t) is given by 

0(t) EE (| E | e -Ht(+fc£+B»?) e -«(-&oJ+Ba?)| | E ) = 

2b 2 

1 ~ sin2 ^ • ( 8 ) 

Here il = \J B 2 + b 2 . p(t) is mixed when \0(t)\ < 1. 
The entanglement with a single spin is reversible because 



there are periodic revivals at the times when 0{t) = 1. 
We note that for b 2 <C B 2 the overlap 0(t) remains 1 
and the qubit state stays pure. 

Even just one environmental spin can lead to irre- 
versible full decoherence when the spin interacts with its 
thermal bath. In the absence of the qubit, thcrmalisa- 
tion of the spin density matrix pe towards |/ can be 
described by a master equation 

T r —pE = <?-Peo-+ + <7 + pe<J- - Pe , (9) 
at 

where cr^ = (cr^ T i<J^ )/2 are spin-1/2 lowering/raising 
operators and T r is the spin relaxation time of the order 
of minutes. This master equation can be formally unrav- 
elled by a stochastic evolution of the spin state. The spin 
state makes a down(up) flip within time dt with a proba- 
bility given by dt/T r times the probability that the spin 
is in the | Te)(| Te)) state. The master equation describes 
evolution of an average over such stochastic trajectories 
of the spin state. 

Let us apply this stochastic unravelling to the state 
of the spin and the qubit. The \ip(t)} in Eq.(^) can be 
rewritten as 

m)) = (l? sinm ) H T>_/31 ^I Te > + 

(rasVm + ^sinQtJ (a| t)+/3| |» | 4e) (10) 

The up transition at the time t can be described as a 
projection of \ip(t)) on (jg | followed by a spin flip 
I 4e) - * I Te)- As a result \ip(t)} jumps to 

|^ + )) = (a|T)+m>)|T E ). (11) 

The resulting state of the system is the same as its initial 
state: this up transition causes no decoherence. 

The probability of the down transition in the interval 
(t, t + dt) is dt/T r times the probability that \i/j(t)) is in 
the | Te) state which is |(Te H(t))\ 2 ~ (b 2 /B 2 ) sin 2 Bt w 
b 2 /2B 2 , where we assume b 2 -^i B 2 . The probability is 
dtb 2 /2B 2 T r . The down transition is a projection on (Te I 
followed by a spin flip | Te) — ► | Ie)- \ip(t)) jumps to 

|^ + )) = HT)s-/3|I)s)II)e. (12) 

The resulting system state has a flipped sign of /3 as 
compared to the initial system state (||). This sign flip 
also flips the sign of the off-diagonal elements in the 
reduced qubit density matrix, compare Eqs.(^,Q). The 
evolution of p(t) including spin relaxation is obtained 
after an average over stochastic trajectories. The off- 
diagonal elements of p(t) average to zero on a timescale 

of Tg = 2T r B 2 /b 2 . This is the time of full decoherence 
due to spin relaxation. We note that for B 2 ^> b 2 the 
time is much longer than T r . 



2 



III. N-SPIN ENVIRONMENT 

Decoherence becomes more efficient when there are 
many copies of the environmental spin numbered by 
i = 1,...,N. Interaction Hamiltonians are iJqj-Ej = 
biOzO^ 1 . The initial product state 
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|V(0)) = (a| T)+/?| W 



(13) 



evolves into a \ip(t)) such that the reduced state of the 
qubit has the form of Eq.(M) but with 



N 



2b 2 



o(t) = J](i-^ S in 2 aa 



(14) 



where Hi — \JB 2 + b 2 . The early time overlap is 0(t) = 
l-2Nb 2 t 2 , where b 2 = £\ h V N - decays from 1 on a 
timescale of = 1/^/Nb. For b\ <C B 2 different factors 
in Eq.jl^) oscillate with different frequencies dispersed 
in the range B ± b 2 /B. For t » 2nB/b 2 the phases of 
oscillating terms in different factors become effectively 
random and the overlap averages to 



O = 1 



Nb 2 



1 



1 

Q 



(15) 



For Q — B 2 /Nb 2 3> 1 the system state remains pure; 
there is only a negligible partial decoherence. Fig. 2 shows 
0{t) for a sample of tVs at two different values of B. 




FIG. 2. The overlap 0(t) from Eq.@ for TV = 10 3 spins 



with b 



1 at two different values of B. bi's were chosen 



at random with a uniform probability distribution between 
and \/3. The spin relaxation is not included in this plot or, 
in other words, T r is assumed to be inifinite. 



Let us include spin relaxation. To get full decoherence 
it is enough that only one out of N spins makes the down 
transition so the decoherence time is N times shorter 
than for a single spin, = (2T r B 2 /b 2 )/N = 2QT r . 
Q 1 makes this decoherence time is much longer than 

T r . 

It is worthwhile to compare the large Q limit with the 
case of B = when Eq.(n4h gives 



JV 



0(t) = Y[cos2bit 



(16) 



The different factors, which oscillate with frequencies in 
a range proportional to b, go out of phase after a time of 
2ir jb when the overlap averages to zero. For B = the 
time t^ 1 ) = 1/y/Nb when 0{t) decays away from 1 is a 
time of full decoherence. 



IV. THERM ALIZED ENVIRONMENT 

So far we considered the single spin environment and 
the multispin environment with a fully polarized initial 
state. A thermal initial state is a weighted average over 
different partially polarized states, like e.g. | Tei)| Ie 2 
) . . . | Te n )- The overlap ( |l4| ) is the same for any such 
partially polarized state as for the fully polarized state so 
the estimates of decoherence time for the polarized state 
go through unchanged for the thermal density matrix. 

The mixedness of the initial environmental state 
does matter when we finally take into account the z- 
component of the qubit magnetic field at the location of 
spins which gives an extra term in the interaction Hamil- 
tonian 



H, 



Q-E 



JV 

E^ 1 



(17) 



This Hamiltonian, when applied to a partially polarized 

spins state, gives -Hq-e = a z J2iLi \^l\Ph where p. t = ±1 
is a polarization of the z-th spin. If p^s were static, then 
this Hq_ e could be measured at the beginning of the 
quantum computation and balanced by a ovcounterterm 
in the qubit Hamiltonian, as sugg ested in Ref. @. The 
problem is that p^s are not static because they flip at 
the rate of one spin per T r /N. The total polarization 
p(t) — Pi(i) performs a random walk with one ±2 

step every T r /N so that [p(t) - p(0)} 2 = 4tN/T r , where 
the over line means here an average over random walks. 
The Hamiltonian Hq_ e w a z b[p(t) — p(0)] is random so 
the qubit states | f) and | |) accumulate opposite ran- 
dom phases ±<^(£), which grow as 4>{t) 2 = ANb 2 ^ /3T r . 
They lead to full dephasing after a dephasing time — 
(S^Tr/Nb 2 ) 1 / 3 . In this discussion we assume Q » 1 so 
that entanglement is negligible. 
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V. NUMERICAL ESTIMATES 

The external magnetic field B is bounded from above 
by the critical field of the type I superconducting material 
H c . This parameter is the highest (0.206 T) for niobium. 
In this respect niobium is much better than alluminium 
which has much lower H c = 0.01 T. We choose niobium 
and set B — 0.1 T, There are materials like lead or 
tin that have abundant isotopes with zero nuclear spin. 
We do not consider them here because niobium and allu- 
minium microtcchnology is more advanced and, what is 
even more important, even lead or tin or their substrates 
would have spin impurities requiring environment engi- 
neering. 

In the qubit design of Ref. |J the 1/im x 1/im super- 
conducting loop is made of an alluminium 0.5/im x 0.5/im 
wire with persistent current of lOOrtA We suggest nio- 
bium instead of alluminium and the loop to be just 0.1/im 
thick. 0.1/xm is roughly twice the penetration depth so 
that B can penetrate through the superconductor. Note 
that the critical field is higher in such thin films than in 
bulk superconductor |1. Our superconducting loop will 
produce a flux of 10 — 4< &o which is an order of magnitude 
less than in the design of Ref. f| . The loop is also thinner 
so it contains fewer spins. The superconductor volume of 
10 _19 m 3 contains 10 10 of niobium nuclear spins. They 
are subject to an average b 2 = 10~ 15 T 2 . The quality 
factor Q = B 2 /Nb 2 = 10 3 gives the overlap O = 0.999. 

If the external field were B = 0, then full decoherence 
would take place after 



-(i) 



1 



Nb 



10" 5 s , 



(18) 



where we assume the nuclear spin magnetic moment of 
10 7 Hz/T. For our B = OAT after we get only neg- 
ligible partial decoherence with O = 0.999. The full de- 
coherence due to spin relaxation could happen after 



.(2) 



2QT r 



10 5 s 



(19) 



where we assume T r to be of the order of minutes, but 
the dephasing time 



.(3) 



/37T 2 T r 



1/3 



10~ 3 s, 



(20) 



is much shorter. Even this 1ms time can be signifi- 
cantly increased by decreasing temperature because T r ~ 
exp(A//cBT), where A is the superconducting gap. Re- 
duction of temperature from 50mK down to 5mK in- 
creases r^ 3 ) 25 times. 



negligible. The same technique can be used for any qubit 
design which involves superposition of states with dif- 
ferent magnetic fields. According to our estimates, the 
present day persistent current qubit proposal || is not 
far from achieving a 1ms dephasing time due to the finite 
spin relaxation time T r . lms is sufficient for the planned 
quantum computation at a GHz rate. To achieve this 
desired decoherence time it is important to replace al- 
luminium by niobium and make the circuit more micro- 
scopic than in Ref. ||. It is a common wisdom that 
more microscopic systems are better from the point of 
view of decoherence, compare e.g. our formula for 
We showed that in addition to this much expected be- 
haviour, there is a threshold below which a sufficiently 
microscopic qubit can be made free from entanglement 
with the spin bath. 
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VI. CONCLUSION 



Maximalization of the quality factor Q = B 2 /Nb 2 so 
that Q > 1 makes entanglement with spin environment 
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